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Abstract 

The dynamic relaxation process for the two dimensional Potts model 
at criticality starting from an initial state with very high temperature and 
arbitrary magnetization is investigated with Monte Carlo methods. The 
results show that there exists universal scaling behaviour even in the short- 
time regime of the dynamic evolution. In order to describe the dependence 
of the scaling behaviour on the initial magnetization, a critical character- 
istic function is introduced. 
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For long it was believed that no universal behaviour would be present in 
the short-time regime of critical dynamics. However, for the critical relaxation 
process starting from an initial state with very high temperature and small mag- 
netization, it was recently argued by Janssen, Schaub and Schmittmann with 
renormalization group methods that there exist universality and scaling even at 
early times, which sets in right after a microscopic time scale tmic- For the 0{N) 
vector model with dynamics of model A, the authors derived a dynamic scaling 
relation which is valid up to the macroscopic short-time regime, 

M^^^ {t, r, mo) = r '^^/'^M^^') {b''t, h^'^r, If^m^) (1) 

where M'^^'^ is fc-th moment of the magnetization, t is the dynamic evolution time, 
r is the reduced temperature, and the parameter h represents the spatial rescaling 
factor. Besides the well known static critical exponents /3, v and the dynamic 
exponent z, a new independent exponent xq which is the anomalous dimension 
of the initial magnetization mo is introduced to describe the dependence of the 
scaling behaviour on the initial conditions. Based on the scaling relation it is pre- 
dicted that at the beginning of the time evolution the magnetization surprisingly 
undergoes a critical initial increase 

M{t)^mot^ (2) 

where the exponent 9 is related to the exponent a;o by 6' = {xq — [3/v)/z. This 
makes the effect of mo very prominent. 

Numerical simulation supports the above predictions for the critical short- 
time dynamics. The exponent 6 for the Ising model was first obtained indirectly 
through the power law decay of the auto-correlation ^ . Recently the initial in- 
crease of the magnetization in was observed for the Ising model and the Potts 
model, and the exponent 9 was directly measured The scaling relation 

and the universality are confirmed [^, |^, |^. The investigation of the universal 
behaviour of the short-time dynamics not only enlarges the fundamental knowl- 
edge on critical phenomena but also, more interestingly, provides possible new 
ways to determine all the static exponents as well as the dynamic exponents from 
the short-time dynamics, either based on the power law behaviour of the observ- 
ables at the beginning of the time evolution H, , or on the finite size scaling 
[ Pm , |TT|. Moreover the universal behaviour of the short-time dynamics is found to 
be quite general, e.g. in the dynamics beyond model A or at the tri-ciritical point 
[0, |T3|, [I^ , in connection to ordering dynamics or damage spreading |T^, ^ and 



on the surface critical phenomena Therefore thorough understanding of the 
universality and scaling for the short-time dynamics is urgent and important. 

The scaling relation ([1|) is valid only under the conditions that the initial state 
is at very high temperature and with small initial magnetization. Are there some 
reasons that the universal behaviour emerge only in the critical relaxation starting 
from such a special initial state? If one believes that the large time correlation 
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length is essential for the universality, the only background one would find is that 
the initial temperature Tq = oo and the initial magnetization mo = are the 
fixed points under the renormalization group transformation. Therefore one may 
wonder whether there exists universal behaviour in the critical relaxation process 
starting from an initial state with very high temperature but initial magnetization 
mo — 1, since mo = 1 also corresponds to a fixed point. Actually the critical 
relaxation process with mo = 1 for the Ising model and the Potts model have 
been investigated with Monte Carlo simulations |T^, |TP|, |n|, |^ . The results show 



that universality and scaling appear to be valid also in the early stage of the time 
evolution. 

In this letter we are more ambitious. We study whether there exists univer- 
sal scaling behaviour in the short-time regime of the critical relaxation from an 
initial state with very high temperature and arbitrary magnetization. If the large 
time correlation length plays an essential role for the universality in the critical 
dynamics as it was pointed out above, the presentation of the universal behaviour 
should not rely on from what initial conditions the critical relaxation starts, and 
only the scaling behaviour of the initial conditions should be considered very care- 
fully. Since the arbitrarily valued initial magnetization is no more around a fixed 
point, a critical exponent Xq is not sufficient to describe the scaling behaviour of 
the initial magnetization. In the language of the renormalzation group method, 
the effective dimension of the initial magnetization will in general depend of the 
initial magnetization itself. In order to describe this phenomena, we introduce 
a critical characteristic function rather than a critical exponent. For the k-th 
moment of the magnetization, the generalized scaling relation may be written as 

M^''\t, T, L, mo) = h-^^/^M^^\h'H, b^/'^r, b'^L, x{b, mo)) (3) 

For convenience of later discussions, finite systems have been considered here 
and L is the lattice size. The scaling behaviour of the initial magnetization mo 
is specified by the critical characteristic function mo), which in the limit 
mo — i> tends to the simple form 6^° mo, but is in general different. Such a 
generalized scaling form is in a similar spirit as that in the correction to the 
scaling, where non-linear effects of an off-fixed point are considered. Here only 
our initial magnetization is a relevant operator and can be far away from the 
fixed point. The ansatz that the exponents (3, v and z do not depend on the 
initial magnetization mo is based on the assumption that the initial conditions 
should not enter the renormalization of the critical system in equilibrium and near 
equlibrium. In other words, if there is a scaling form in the short-time regime of 
the dynamic relaxation process, the scaling form should smoothly cross over to 
that in the long-time regime. This greatly simplifies the short-time behaviour of 
the critical dynamics. In the neighbourhood of mo = 1 the critical characteristic 
function mo) may also be characterized by a critical exponent. One can also 
realize that in the limit of 6 = 0, x^i'^o) ~^ O5 ^ind when b approaches infinity, 
x{b,mo) — > 1. In this letter we are interested in the more general case, i.e. mo is 
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between and 1, and 6 is a reasonable finite number. The limiting cases will be 
discussed in detail elsewhere. 

Here we stress that the scaling relation in is not trivial even though the 
critical characteristic function mo) looks not so "simple" as a critical ex- 
ponent. The scaling relation relates the time evolution of the observables with 
different initial magnetizations to each other and represents the self- similarity 
of the dynamic systems. All the physical observables as functions of t, r and L 
should be described by the same critical characteristic function mg). Besides 
this, the physical observables are universal functions of the variables t, t and L 
up to a non-universal scaling constant and the dependence of the observables as 
well as "^o) on mo is expected to be universal up to a rescaling of the initial 
magnetization mo. 

As a concrete example, we consider the two dimensional Potts model, for 
which a quite accurate dynamic exponent z has been obtained from the universal 
behaviour of the short-time dynamics 0. The Hamiltonian for the q state 
Potts model is given by 

H = KY.6^^^,^, ai = l,...,q (4) 

<ij> 

where < ij > represents nearest neighbors. In our notation the inverse tem- 
perature is already absorbed into the coupling K. It is known that the critical 
points locate at Kc = log{l + ^/q). In this paper the three state [q = 3) case 
will be investigated. In generating the random initial configurations, we have 
sharply prepared the initial magnetization in order to avoid extra finite size ef- 
fects from the fluctuation of mo in the finite systems. After the preparation of the 
initial configurations the system is released to evolve according to the Heat-bath 
algorithm . 





mo 


0.14 


0.40 








Ud 




Ud 


(36, 72) 
(72, 144) 


x(2,mo) 


0.1807(05) 
0.1800(04) 


0.1804(08) 
0.1798(04) 


0.5298(06) 
0.5307(06) 


0.5300(36) 
0.5302(45) 


(36, 144) 


x(4,mo) 


0.2328(09) 


0.2324(11) 


0.6918(28) 


0.6910(78) 



Table 1: x(6, mo) measured from the magnetization and the Binder-type cumu- 
lant. 

For simplicity, we take r = and therefore the exponent u will not enter the 
calculation. The exact value of the static exponents (3/i> = 2/15 and the dynamic 
exponent z = 2.196 obtained from the power law decay of the auto-correlations 
[05 01 ^ill be taken as input. To verify the scaling relation (0) and determine the 
critical characteristic function x(6, mo), we perform the simulation for a pair of 
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lattice sizes Li and L2 and measure the time evolution of the magnetization and 
the second moment defined as 

MW(t,mo) =< Ei^a.m - I)]' >, k = l,2 (5) 

i 

where the average is taken over independent initial configurations and the ran- 
dom forces. In order to reduce extra errors from M^''\t, 1), especially when mg 
becomes bigger, we introduce a magnetization difference 

Ma{t, mo) = M(i) (t, 1) - M^^) {t, mo) (6) 

and a Binder-type cumulant 

From the scaling coUaps of Mdit^rriQ) or Ud(t,mo) for two lattices with suitable 
initial magnetizations, we can estimate the values of the function mo) at 
b = L2/L1 for different mo- 

In Fig. 1 the scaling plot for the magnetization is displayed for the lattice 
L2 = 144 with initial magnetization mo2 = 0.14 and the lattice Li = 72 with 
suitably correponding initial magnetizations moi. The stars represent the time 
evolution of the magnetization of the lattice L2 = 144. The crosses are the same 
data but rescaled in time and multiplied by an overall factor 2^^'^. If one can 
find a moi for lattice Li = 72 such that its time evolution of the magnetization 
fits to the crosses, the scaling is valid and from the scaling relation (|^) one gets 
x(2, mo2) = Tnoi. Practically we have performed the simulation for Li = 72 with 
two different initial magnetizations mo = 0.175 and mo = 0.185, for which the 
magnetizations have been plotted by the lower and upper solid lines in Fig. 1. 
By linear extrapolation, we obtain the time evolution of the magnetization with 
initial values between these two values. Then we can estimate moi. The solid 
line laying on the crosses in Fig. 1 is the time evolution of the magnetization for 
L = 72 with an initial magnetization mg = 0.1800(4) which has the best fit to the 
crosses, i.e. x(2,0.14) = 0.1800(4). From the previous numerical simulations for 
the short-time dynamics [|, ||, 0], we know that the microscopic time scale tmic 
is negligibly small for the Heat-bath algorithm, at least for the measurement of 
the dynamic exponent 6 or the exponent Xq. In our calculation we have carried 
out the fitting procedure in a time interval of [10,200] in the time scale of lattice 
L = 72. From Fig. 1 one can also see clearly that the scaling relation is valid 
starting from the very early stage of the time evolution. 

In Fig. 2 a scaling plot is shown for the lattice L2 = 144 with mo2 = 0.14 
and Li = 36. From a fitting in a time interval of [10,80] we obtained x(4, 0.14) = 
0.2328(9). In a similar way, x{b, ''^0) may independetly obtained from the scaling 
coUaps of the Binder-type cumulant Uci(t,mo). Altogether we have performed 
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Figure 1: The scaling plot for the magnetization with (Li, L2) — (72, 144). 



the simulation for mo2 = 0.14, 0.22, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80 and 1.00. 
In Tab. 1 the results for two typical values mo2 = 0.14 and 0.40 are given. 
We can see that mo) estimated from both Md{t,mQ) and Ud{t,mQ) are very 
well consistent. In order to see the finite size effect, we have also performed the 
calculation for other pairs of lattices as e.g. (36, 72). In Tab. 1 we can see that 
the finite size effect for the lattice pair (Li, L2) = (36, 72) is already quite small. 
In the simulations, the statistics of L = 36, L = 72 and L = 144 are respectively 
80 000, 40 000 and 8 000. In Tab. 1, errors are estimated by dividing the data 
into four groups. In principle the measurements from the magnetization are more 
reliable than those from the Binder-type cumulants where higher moments are 
involved, at least when the static exponent jS/u is known. The magnetization is 
self- averaging, but the Binder-type cumulant is not. When tuq is getting bigger, 
Udit.mo) is more fluctuating. 

In order to get a more direct understanding of the full critical characteristic 
function x(6, mo), we define an effective dimension x{b, mo) of the initial magne- 
tization mo by xi^^iT^o) — b ^(''■"^o) mo- By the definition naturally x{b, 0) = xq. 



5 



0.70 




Figure 2: The scaling plot for the magnetization with (Li,L2) = (36, 144). 

From the values of x{b, rno), we can calculate x{b, mo). Taking the results from the 
magnetization difference with the lattice pair (Li,L2) = (36, 144) and (72, 144), 
the corresponding effective dimensions x{b, mo) are plotted in Fig. 3. It clearly 
shows that x{b, mo) is a non-trivial function. For example, the values x{2, 0.14) = 
0.363(3) and x(2, 0.40) = 0.408(2) apparently differ from xo = 0.298(6) §. When 
mo varies from zero to one, the effective dimension x(6, mo) first increases and 
then decreases to zero. However, when mo (or x{b,rno)) is approaching mo(or 
x(6, mo))=l it is not the best choice to determine x{b,mo) directly in the way 
discussed in this letter since the dependence of the physical observables on mo 
becomes weaker and we face big statistical fluctuations. 

Finally let us have some more understanding of the dynamic system discussed 
above. Supposing mo is also a 'coupling' of the system, in the ^K—mQ plane of the 
couplings there exists a critical line K = Kc 'ia the sense that the time correlation 
is divergent. In the neighbourhood of this critical line, the time correlation length 
depends only on the coupling K, and its scaling behaviour is characterized by 
the exponent vz. If the critical line is a line of the fixed points, the exponent 
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Figure 3: The effective dimension obtained from the magnetization with 
(Li,L2) = (72,144) and (36,144). Circles are of x(2, mo) and stars are of x(4, mo). 
The hnes are drawn to guide the eyes. 

V together with z and /3 is sufficient to describe the critical dynamic system. 
However, in our case the critical line is not a line of the fixed points. Therefore 
a critical characteristic function should in general be introduced to describe the 
scaling behaviour. It is actually interesting to sec what happens for a system 
in the equilibrium where a critical line exists but the line is not a line of fixed 
points. 

In conclusion, we have numerically simulated the universal short-time be- 
haviour of the dynamic relaxation process for the two dimensional critical Potts 
model starting from an initial state with very high temperature and arbitrary 
magnetization. The results show that the traditional scaling relation should be 
generalized. A critical characteristic function is introduced to describe the scal- 
ing behaviour of the initial magnetization. We demonstrate how to determine 
numerically the critical characteristic function. The study of the short-time dy- 
namics is not only conceptually interesting but also pratically important since it 
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is possible to obtain the static exponents /?, v and the dynamic exponent z of 
the critical systems far before the dynamic process reaches the equilibrium. It is 
challenging to derive analytically the generalized scaling relation in as well 
as the critical characteristic function mo). The application to the dynamic 
field theory, e.g. the stochastic quantization of the field theory, is attractive since 
such a knowledge would be important for the numerical simulation of the lattice 
gauge theory. 
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